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We discuss fermionic zero modes in the two-dimensional chiral p-wave superconductors. We show 
quite generally, that without fine-tuning, in a macroscopic sample there is only one or zero of such 
Majorana-fermion modes depending only on whether the total vorticity of the order parameter is 
odd or even, respectively. As a special case of this, we find explicitly the one zero mode localized on 
a single odd- vorticity vortex, and show that, in contrast, zero modes are absent for an even- vorticity 
vortex. One zero mode per odd vortex persists, within an exponential accuracy, for a collection of 
well-separated vortices, shifting to finite ±E energies as two odd vortices approach. These results 
should be useful for the demonstration of the non-Abelian statistics that such zero-mode vortices 
are expected to exhibit, and for their possible application in quantum computation. 



PACS numbers: 74.20.Rp, 05.30.Pr, 03.67.Lx 

Recently [1-3] there has been considerable interest in 
the structure of fermionic zero modes localized on vor- 
tices of a chiral spinless two-dimensional superconductor 
characterized by p x + ip y order parameter. In part, it is 
stimulated by a proposal [1] that a ground state of such a 
superconductor (for a positive chemical potential) is sim- 
ilar to the Moore- Read (Pfaffian) quantum Hall state [4] , 
thought to describe the v = 5/2 quantum Hall plateau 
[5]. Vortices (corresponding to the Laughlin quasihole- 
like excitations [6] in the Moore-Read state) in such a 
superconductor are thus expected to exhibit a degener- 
ate set of zero modes separated from all other states by 
a gap, and to obey non-Abelian statistics [7], that may 
make them useful for a realization of a "topological quan- 
tum computer" [8] free of decoherence. 

Many of the properties of these zero modes for a single 
vortex have already been discussed in the literature [1- 
3]. However, in our view an explicit discussion of the 
fate and robustness of the zero modes to, for example, a 
local deformation of the order parameter or in the pres- 
ence of many vortices has not appeared in the literature. 
Such questions are of particular interest in view of recent 
proposals for experimental realization and manipulation 
of such non-Abelian states in two-dimensional supercon- 
ductors, such as Sr2RuC>4 [9], the v = | plateau of the 
quantum Hall effect [10-12], and p-wave resonantly in- 
teracting atomic superfluids [2, 13]. 

In this Letter we show quite generally, that for a macro- 
scopic sample (i.e., ignoring the boundary physics), with- 
out fine-tuning, strictly speaking there is only one or zero 
Majorana-fermion mode depending only on whether the 
total vorticity of the order parameter (in elementary vor- 
tex units of 2-7r) is odd or even, respectively. For a col- 
lection of well-separated vortices, within an exponential 
accuracy one zero mode per an isolated odd- vorticity vor- 
tex persists. As two of such vortices are brought closer 
together the corresponding pair of "zero" modes splits 
away to finite ±E (vortex-separation dependent) ener- 
gies. Generically, even- vorticity vortices do not carry any 



zero modes. 

Even in the odd-vorticity case, zero modes only exist 
for a positive chemical potential fi > 0, consistent with 
the existence (absence) of a topological order in a weakly- 
( strongly-) paired ground state of a p-wave superconduc- 
tor stable only for fj, > (jt < 0) [1, 14]. While a p-wave 
superconductor in a solid state context naturally obeys 
fi > 0, in a Feshbach resonant atomic p-wave superfluid 
a chemical potential can be adjusted to be positive via 
an external magnetic field [13], a "knob" that can also 
be used to drive a topological quantum phase transition 
between a strongly- and weakly-paired superfluid ground 
states [1]. 

As a demonstration of a specific realization of this gen- 
eral connection between parity of vorticity and a number 
of zero modes, we compute the eigenfunction of the one 
zero mode localized on a single isolated odd-vorticity vor- 
tex, and show that zero modes are absent for an even- 
vorticity vortex. This symmetric vortex result is in agree- 
ment with a recent study in Ref . [3] , but does not rely on 
a linearization of the fermion dispersion around a Fermi 
surface, and thereby allows us to access the nondegener- 
ate (low chemical potential) regime realizable in tunable 
(via a Feshbach resonance [13]) atomic gas experiments. 
Our results then imply that such zero-modes, residing 
on isolated elementary vortices are always shifted to fi- 
nite ±E energies when an even number of them come 
into proximity [1], with possible deleterious implications 
for a realization of non-Abelian statistics and quantum 
computation. 

To demonstrate these results we begin by first 
discussing the properties of the solutions of generic 
Bogoliubov-de-Gennes (BdG) equations arising in a con- 
text of any superconductor. These coupled Schrodinger 
equations follow from the following Bardeen-Cooper- 
Schrieffer (BCS) Hamiltonian 

H = (">-.."., ~ ajhijol + aiAijaj + a} Ay at) , (1) 
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where indices i, j label space (and in a spinful case, spin) 
coordinates of the fermion creation and annihilation op- 
erators a t , dj . Their canonical anticommutation relations 
ensure that is an antisymmetric operator. Since H 
must be hermitian, so is hij, and the problem is equiv- 
alent to a study of the spectrum and eigenstates of a 
matrix 
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This matrix possesses the following important symmetry 
property 



o\H.o\ = —TL*, 



(3) 



where a\ is the first Pauli matrix acting in the 2 by 2 
space of the matrix TL, Eq. (2). In the terminology of 
Ref. [15], the matrix TL is said to belong to the symmetry 
class D. As a result of this property, it can be seen from 



Wo-iV* = -<riH*ip* = -Eanp* 
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that if tp is an eigenvector of such TL with the eigenvalue 
E, then aiip* is guaranteed to be an eigenvector with the 
eigenvalue —E. As a result, all nonzero eigenvalues of TL 
come in ±E pairs. A special role is played by the zero 
eigenvectors of this matrix, referred to as zero modes. If 
tp is a zero mode, a\tp* is also a zero mode. Taking linear 
combinations tp + u\tp* , i(tp — a\tp*) of these degenerate 
modes, we can always ensure the relation 



a\tp* = tp 



(5) 



for every zero mode. In the absence of other symme- 
tries of TL it is quite clear that generically there is noth- 
ing that protects the total number N z of its zero modes 
under smooth changes of the Hamiltonian matrix that 
preserve its BdG form. However, since non-zero modes 
have to always appear in ±E pairs, as long as the sym- 
metry property (3) is preserved by the perturbation the 
number of zero modes can only change by multiples of 2. 
Thus, while the number N z of zero modes of the Hamil- 
tonian (2) may change, this number will always remain 
either odd or even, with (— 1)^* a "topological invariant" 
[16, 17]. 

The value of this invariant is easy to establish if one 
observes that TL is an even sized matrix, with an even 
number of eigenvalues. Since the number of non-zero 
modes must be even, this implies, quite generally that 
the number of zero modes is also even, ( — 1) N * = 0, and 
strictly-speaking the BdG Hamiltonian docs not have any 
topologically protected zero modes. Furthermore, since 
zero modes must appear in pairs, there can only be an 
even number of accidental zero modes, which will never- 
theless be generally destroyed by any perturbation of TL 
(preserving its BdG structure Eq. (2)). We believe this 
observation was first made by N. Read [17]. 



The situation should be contrasted with that of the 
Dirac operators T>. Those operators, being generally of 
one of the chiral classes in the terminology of Rcf. [15], 
obey the symmetry 

a 3 Va 3 = -V. 

Thus if tp is an eigenvector of T> with the eigenvalue E, 
a 3 tp is an eigenvector with the eigenvalue — E. Thus, 
(after a suitable diagonalization) the zero modes of V 
must obey the relation 
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Namely, they are eigenstates of the a 3 operator, with 
the "left" zero modes ipL coming with the eigenvalue +1, 
and the "right" zero modes tpn labelled by the eigenvalue 
— 1. As the operator V is deformed, the number of zero 
modes changes, but the non-zero modes always appear 
in pairs, where one of the members of a pair has to be 
"left" and the other "right" . Therefore, while the number 
of zero modes is not an invariant, the difference between 
the number of left and right zero modes is a topological 
invariant, determined (through the index theorem) by 
the monopole charge of the background gauge-field. 

Contrast this with zero modes of TL, which obey the re- 
lation Eq. (5). Because of the complex conjugation on tp, 
these zero modes cannot be split into "left" and "right" . 
Indeed, even if we tried to impose a\tp* — —tp, a simple 
redefinition of tp — > itp brings this relation back to Eq. (5). 
Thus, the most an "index theorem" could demonstrate 
in the case of the BdG problem, is whether there is or 
exactly 1 zero mode. Moreover, since the BdG problem 
is defined by an even-dimensional Hamiltonian, generi- 
cally there will not be any topologically protected zero 
modes [16, 17]. 

Yet it is quite remarkable that in the case of an iso- 
lated vortex of odd vorticity in a macroscopic sample 
(i.e., ignoring the boundaries) of a p x + ip y superconduc- 
tor of spinless fermions, there is exactly one zero mode 
localized on this vortex [1-3, 18]. To be consistent with 
above general property of the BdG Hamiltonian (namely, 
that the total number of zero modes must be even) an- 
other vortex is situated at the boundary of the system 
[1, 17] , preserving the overall parity of the number of zero 
modes. Hence, although even in this odd-vorticity case 
the one zero mode is not protected topologically, able to 
hybridize with a vortex at a boundary of the sample, it 
survives (up to exponentially small corrections) only by 
virtue of being far away from the boundary (and from 
other odd-vorticity vortices). 

To see this explicitly we consider the BdG equations 
for a two-dimensional p x + ip y superconductor 

d 
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Here A(r) is the order parameter of the superconductor, 
z = x + iy, z = x — iy are the two-dimensional complex 
coordinates, m is the fermion mass, and \i is the chemical 
potential. Eq. (6) is of course a particular case of the 
eigenvalue equation for a matrix of the form given in 
Eq. (2), with the vector tp represented by 



(7) 



For a uniform (vortex-free) order parameter, A(r) = A , 
it is easy to solve Eq. (6) in terms of plane waves, finding 
the spectrum 
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Since has a gap for all k (with the exception of the 
critical point at [i — [1, 14]), consistent with above 
discussion, there are no zero modes of Eq. (6) in the 
absence of vortices. 

Now consider a superconductor with a symmetric vor- 
tex of vorticity I. The order parameter is then given by 
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where r, <p are the polar coordinates centered on the vor- 
tex and f(r) is a real function of r that vanishes at small 
r. Then the BdG equations take the form 
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Next we observe that for the case of a vortex of even 
vorticity, £ = 2n, we can eliminate the phase dependence 
of Eq. (10) entirely . Indeed, making a transformation 
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leads to equations 
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0), and smoothly deforming f(r) into a constant equal 
to its asymptotic value at large r; in order to be smooth, 
the deformation must preserve the BdG structure Eq. (2) 
and the vorticity of the order parameter. These equations 
then become equivalent to Eq. (6) for a constant, vortex- 
free order parameter A(r) = A with an exact spectrum 
Eq. (8), that for fj, ^ clearly does not exhibit any zero 
modes. 

As Eqs. (12) are smoothly deformed, in principle it 
is possible that for a particular deformation some of its 
eigenstates will become zero modes (although, as demon- 
strated above, this can only happen in ±E pairs, leading 
to an even number of these). However, these modes will 
not be topologically protected, and even a small defor- 
mation of, say, the order parameter shape f(r) will de- 
stroy these modes. We note that this argument easily 
accommodates vortices that are not symmetric, as those 
can be smoothly deformed into symmetric ones without 
changing the topologically protected parity of N z . The 
conclusion is that generically there are no zero modes in 
the presence of an isolated vortex of even vorticity. 

In fact, if any doubts remain, it is also possible to di- 
rectly demonstrate the absence of zero modes in Eq. (12), 
simply by following the arguments parallel to those given 
after Eq. (20). However, the arguments presented above 
are more general and robust, and can be used to establish 
the claim even for non-symmetric even-vorticity vortices. 

The situation is drastically different if the vorticity of 
a vortex is odd, i.e., £ = 2n — 1. Indeed, in that case 
the transformation Eq. (11) cannot entirely eliminate 
such vortex from the equations (even with the help of 
a smooth deformation), leaving at least one fundamental 
unit of vorticity. This thereby guarantees at least one 
zero mode localized on the odd-vorticity vortex. To see 
this, recall that due to the condition Eq. (5) together 
with the definition Eq. (7), the zero mode satisfies 
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Combining this with the transformation Eq. (11), we find 
the equation for the zero mode 
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Now we note that these equations are topologically equiv- 
alent to the BdG equations without any vortices. Indeed, 
the only difference between these equations and those 
for a uniform condensate is the presence of the terms 
2in/r 2 [d/d(p], n 2 /r 2 , and /(r) that is a constant at large 
r and vanishes in the core of the vortex for r < r coro . We 
can imagine smoothly deforming these equations to get 
rid of the first two terms (for example, by replacing them 
with a (n 2 jr 2 — 2in/r 2 [d/dp\) u and taking a from 1 to 



We look for a solution to this equation in terms of a 
spherically symmetric real function u(r). This gives 
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A transformation 
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brings this equation to the more familiar form 
x" x' J f 4 (r) n 2 \ 

This is a Schrodinger equation for a particle of mass m 
which moves with angular momentum n in a potential 
mf 4 /8 that is everywhere positive. We observe that this 
potential vanishes at the origin, and quickly reaches its 
asymptotic bulk value /o away from the origin. Then for 
[i > m/o/8, there always exist a solution to this equation 
finite at the origin and at infinity. Moreover, if fi < 
m/o/8, then the solution finite at the origin will diverge 
at infinity as 

x ~ eV^- 2 " 1 ". (18) 

Combining this with Eq. (16), we observe that u(r) re- 
mains a bounded function at infinity as long as /j, > 0. 
Thus the conclusion is, there exist a zero mode as long 
as (i > 0. For the special case of the n = vortex of 
vorticity — 1, the small and large r asymptotics of the 
solution we found here was discussed recently in Ref. [2] . 

In the simplest London approximation of a spatially 
uniform condensate with f(r) — f for all r except inside 
an infinitcsimally small core, the zero mode localized on 
an isolated odd-vorticity vortex is simply given by 



)J n yry - m 2 ^fj e tt/o 2 ^ for fj, > m^-, 
I n (r\Jm 2 & - 2nif?j e~^ / ° r , for < ^ < m^f, 

(19) 

where J n (x), I n {x) are Bessel and modified Bessel func- 
tions. 

We note that it may seem possible to construct ad- 
ditional zero modes in the following way. Instead of 
the ansatz of a rotationally invariant u(r) just before 
Eq. (15), we could have chosen an ansatz 

u(r, <p) = u a (r) e iav + u_ a (r) e~ iav . (20) 

Then two second order differential equations follow re- 
lating these two functions. Generally there are four so- 
lutions to these equations. Boundary conditions at the 
origin r — select a subset of two of these solutions. 
Boundary conditions at infinity select a different subset 
of two solutions. However, barring a coincidence, none of 
those solutions finite at the origin are also finite at infin- 
ity. Even if such coincidence arises for some special value 
of fi, by the above arguments, the additional zero modes 
must appear in topologically unprotected pairs, that will 
be split to finite ±E energies by a slight generic deforma- 
tion of the potential (order parameter distortion). Hence 



we conclude that generically there will be no additional 
zero modes (except the one found above) for an odd- 
vorticity vortex. 

Thus we indeed find that the number of zero modes 
in a symmetric odd-vorticity vortex must be one. Since 
a smooth deformations of the order parameter can only 
change the zero mode number by multiples of two, an 
arbitrarily shaped odd-vorticity vortex must have an odd 
number of zero modes. However, any number of zero 
modes other than one is not generic and will revert to one 
under an arbitrary deformation of the order parameter. 

To summarize, the results presented here establish the 
robustness of the zero modes localized on well-separated 
(''separation 3> 1/(toA)) odd-vorticity vortices and sup- 
port the idea that they can eventually be used to demon- 
strate non-Abelian statistics and perhaps even for quan- 
tum computation. 
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